Quantum Dots with Disorder and Interactions: A Solvable harge-g Limit 



Ganpathy Murthy^ and R. Shankar^ 
^Physics Department, University of Kentucky, Lexington KY 40506-0055 
^ Department of Physics, Yale University, New Haven CT 06520 
(February 1, 2008) 

We show that problem of interacting electrons in a quantum dot with chaotic boundary conditions 
is solvable in the g ^ oo limit, where g is the dimensionless conductance of the dot. The critical 
point of the g — oo theory (whose location and exponent are known exactly) that separates strong 
and weak-coupling phases also controls a wider fan-shaped region in the coupling-l/(? plane, just as 
a quantum critical point controls the fan in at T > 0. The weak-coupling phase is governed by the 
Universal Hamiltonian and the strong-coupling phase is a disordered version of the Pomeranchuk 
transition in a clean Fermi liquid. Predictions are made in the various regimes for the Coulomb 
Blockade peak spacing distributions and Fock-space delocalization (reflected in the quasiparticle 
width and ground state wavefunction) . 
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The union of interactions and disorder in electronic 
systems poses a nasty problem: Techniques that work 
when only one or the other is present fail when they 
coexist. Here we discuss a class of problems involving 
quantum dots (QD's) where the evil twins can be tamed, 
thanks to a small parameter 1/g, g being the dimension- 
less conductance of an open QD. The problem is inter- 
esting for experiments in the Coulomb Blockade (CB) 
regime^"'*. 

We consider d ~ 2 dots characterized by Ep, the Fermi 
energy, Et , the Thouless energy (which is the amount by 
which the uncertainty principle broadens the electronic 
energy levels in the time it takes to traverse the dot) , and 
A, the average single particle level spacing, with Ep >> 
Et >> a. Et also measures the band around the Fermi 
energy wherein the energy levels and wave functions {g 
in number) may be described statistically^ by Random 
Matrix Theory (RMT)^'^. The randomness here is due 
to the chaotic boundary conditions, with the mean free 
path I equal to the sample size L, and Et — Et ~ hvp/L, 
where Er is the level width due to scattering. 

Our hamiltonian is: 

^ = Y1 ^i^o^ea + \Y1 Vc.fijsi^ii'li^j'fps (1) 

where Eq are single-particle levels that obey RMT statis- 
tics, have a mean spacing A and range from — gA/2 to 
gA./2. In the following we will supress spin for simplicity, 
pointing out how its restoration modifies various results. 

Choices for Vafs-yS range in previous work from all of 
them being independent gaussian variables* to the Uni- 
versal Hamiltonian^"^ ^, wherein Vappa = uq A, (all oth- 
ers zero). The sole parameter uq clearly couples to Q^, 
Q being the total charge. For the spinful case a term 
coupling to total is also included. 

We employ the choice made by Murthy and Mathur^^ 
(MM) who appeal to the Renormalization Group (RG) 
approach developed by one of us^'^, where one integrates 
out modes that lie outside a narrow band centered on 



the Fermi energy to expose the low-energy physics. In a 
clean system the RG leads to Landau's Fermi Liquid (FL) 
parameters^'* as fixed-point couplings^^. With disorder 
one expects this RG to work till we come down to E^. ~ 
Et, still leading to FL couplings since Et << Ep. At 
this point disorder is included exactly by switching from 
the k basis to the disorder basis: 

A 

J2u{e - 0') [C(k)0^(k') - C(k')'^Mk)] * - '^7) 

kk' 

where 9 and 6' are the angles of the momentum vectors 
k and k' and Vap-ys is just the transcription of the Fermi 
liquid interaction VpL — -f^ Skk' "(^ ~ ^'O^k'^k'- We 
will Fourier resolve the interaction as u(9 — 9') = uq + 
J2m=i ''^m cos [m{9 — 6')]. The constant uq controls V's 
with nonzero average {Vaf3f3a) = "(^: = "oA. 

All couplings have nonzero fluctuations given by 

A^ 

{^af3js) - (Vq/375)^ = ^ ^ "m- (2) 
^ m=l 

Couplings with zero average are discarded in the Uni- 
versal Hamiltonian^"^*. Can these couplings be relevant 
in the RG sense and dominate low-energy physics, de- 
spite their small size in the original Hamiltonian? To 
answer this, MM*^ resorted to further fermionic RG*^ 
within Et- Their one- loop /3- function yielded a critical 
point at it^ = — 1/ ln2 for m > 0. (Upon including spin 
the instabilities can occur in the charge or spin channels, 
the critical coupling being —1/2 In 2 for each channel.) 
For u„i > — l/ln2, the flow led to the Universal Hamil- 
tonian, while for Um < — l/ln2 there was a runaway to 
strong coupling. The two phases seemed clearly related 
to two types of behavior in level spacings and the nature 
of quasiparticles. However, uncertainty surrounding the 
nature of the strong-coupling phase and even its very ex- 
istence, predicated as it was on a one-loop calculation and 
a fixed point of order unity, impeded further progress. 
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FIG. 1. Diagrammatic proof of the large-A'^ (with N = g) 
nature of the theory. Disorder-averaging internal hnes forces 
1 = m with a free sum. 

We have verified that the theory is solvable in the limit 
g — !■ oo, with l/g playing the role ofl/N in large-N theo- 
ries. The key idea is illustrated by proving that the four 
point function Gafi-yS is just a sum of repeated particle- 
hole bubble diagrams. Recall that in theories with inter- 
action V — A/iV(^. V'j'/'j); only such iterated 
bubbles survive, since only they have a free sum over N 
for each extra loop. The situation here is similar. Con- 
sider the one-loop diagrams in our theory. Each vertex 
carries a sum over two momenta ( k and k' in Fig. la 
for the bare vertex), and each propagator is diagonal in 
the disorder eigenvalue index, but not the k index. In 
Fig. lb, the internal states involve a sum over terms 
of the form (/)* (l)</)^(m)0* (m)(/)^(k) Replacing the sum 
(over /i, I/, 1, m ) by the ensemble averages using 



(0;(ki)(/)^(k2)C(k3)'^.(k4)> 



6126: 
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0(1/5^) (3) 



(and neglecting fluctuations down by l/g) we find that 
m = 1 and that there is a free sum over g values of either 
one of them. (An exception occurs when either vertex 
involves uq which does not flow^^.) The reader may check 
that other possible one-loop particle-hole and particle- 
particle diagrams are down by l/g because one or more 
of the external k labels creep into the diagram. The same 
logic holds for higher orders. Consequently, MM's one- 
loop /3- function and critical exponent v = l/(3'{u*) = 1 
are exact and the phase transition is real. We now ask 
what it corresponds to. 

Let us consider just one Um < and factorize: 
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— \um \ cos{m9 ~ m9') 
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dcr exp 



2\u^\A 

where cr = (cr^ , cr^) has two components, as does M: 
(M\m2) = ^(/)* (k)0^(k) (cos me* , sin m6l) (4) 

k 

Integrating out the fermions we get an effective action 
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-Tr\n[{iuj~ea)I~<T-M] (5) 



where / is the unit matrix. The quadratic part of the 
action is 
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upon using the disorder-averaged result 
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where Na is the Fermi occupation factor of level a. The 
&^ term is valid for frequencies well below A. (We shall 
not be interested in faster motion of the collective mode.) 

To reconfirm the \a,Tge-g nature of the theory, one de- 
fines & = cr/g and evaluates the Tr In, term by term, 
and finds that S has a g^ in front, which plays the role of 
1/h. If 1/g'^ =" h" = 0, we are in the classical limit and 
spontaneous symmetry breaking is possible (even) for a 
single degree of freedom, while if 1/g^ > 0, this is impos- 
sible due to "quantum" fluctuations. If we include more 
Mm's, the corresponding cr's will make additive contribu- 
tions to S. As soon as one of them breaks symmetry, the 
rest will not matter. 

The phase diagram is shown in Figure 2. Let us first 
focus on l/g = 0. For m,„ > u^, (cr) — 0. For u,n < 
M^, (<t) ^ with a Mexican hat minimum located by 
balancing the quadratic term we have above with the rest 
of S. This leads to {\cr\) = Cg\/]r[ where C is a number 
of order unity and r — u^^^ — (u^)~^. At l/g = there 
are no fluctuations and cr can sit anywhere on a circle. 

What happens at l/g > follows from its role as the 
pre-factor in the action of a: the critical point can be 
felt within a critical "fan" just as a T = critical point 
can influence T > physics^^. The fan is F-shaped since 
both r and l/g scale with exponent 1. 

For l/g > the symmetry of the Mexican hat valley 
is explicitly broken by sample-specific corrections to next 
order in l/g yielding a unique minimum at some angle 
9. The stability of this minimum to "quantum" fluctu- 
ations (at nonzero l/g) is examined by reading off the 
hamiltonian for cr from Eqn.(6) for the effective action. 
The angular part of cr is governed by 
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where L is the angular momentum conjugate to 9, and 
factors other than g are suppressed. The radial part of 
< cr >~ g which appears in both kinetic and potential 
terms is represented by its average value of order g. The 
potential term dominates at large g and ensures that the 
wavefunction is localized near the minimum of V. Unlike 
in the two phases, the critical point (and fan region) has 
very large fluctuations in cr: From the hamiltonian 



= gP' + — 



(9) 



we estimate the ground state wave function to have a 
width 0{g^^^) in a. As the interaction increases we cross 
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FIG. 2. Schematic showing the critical point at 
— — l/log2 in the g oo limit on the x-axis. As g 
decreases one moves in the y-direction and the critical point 
controls the behavior of all physical quantities in a "fan" . 

over from the Universal Hamiltonian phase through the 
critical region to the symmetry-broken phase. 

What does the order parameter correspond to? In a 
pure system when any Landau parameter Um < — 2 ( 
< — 1 with spin), the Fermi surface undergoes a Pomer- 
anchuk shape transition^'*' For example if m = 2 the 
shape is an oval*^, with the director being a Goldstone 
mode. The present transition is the disordered version 
of this, with small disorder terms making the Goldstone 
mode massive. (An anisotropic gap at the Fermi sur- 
face could also arise when considering instabilities in 
the spin channel^^). Note that the phenomenology de- 
scribed above can occur on top of a mesoscopic Stoner 
transition^^'^'*°. 

The phase diagram clarifies at least two issues. First, 
we can relate the fluctuations of cr to "delocalization 
in Fock space" which concerned the crossover of the 
quasiparticle width to a Breit-Wigner form as a function 
of its energy. Note that a acts as a background field for 
the fermions in the Hartree -Fock Hamiltonian 



HF 



(10) 



Outside the critical fan Hh f is controlled by one sharply 
peaked value of < cr > (either zero or somewhere on 
the Mexican hat's valley) and states are descendants of 
a single Slater determinant. In the fan cr suffers large 
critical fluctuations. Thus we expect the ground state 
to go from being localized to delocalized and back to 
localized as we raise the coupling at fixed g. This re- 
entrance has been seen in numerics^*. We connect the 
quasiparticle width^° to its energy e, by arguing that e 
plays the same role as 1/g or T, and that there will be 
a fan in the u — e plane. As the quasiparticle energy 
increases in either phase, one hits the crossover fan with 
its delocalized states. 

Second, we can understand the dependence of the CB 
peak spacing distribution on m^^'^^. Hhf (Eqn.(lO) is a 
sum of two random Hamiltonians whose widths can be 
added in quadrature to yield an effective level spacing 

A'= VA2+6<(T>2 /g2 (11) 

In the weak-coupling phase < cr >= and there is no 
memory of In the broken symmetric phase we have 



< cr >= CgAy/r, and the effective level spacing increases 
with r, and thus Um- There is a crossover between these 
two regimes in the critical fan. 

Let us now briefly analyze the relevant experiments. 
The Sivan et al^ and Patel et aP experiments are done 
on gated GaAs 2DEG samples, for which r<; w 1 — 1.2 in 
the bulk. Using the area of the sample and the fact that 
these dots are in the ballistic limit we can find both A 
and Et, and thence find 5 w 7 — 14. Sivan et al^ find that 
the CB peak spacing is about 5 times broader than that 
predicted by the Universal Hamiltonian, which describes 
the weak coupling region of our phase diagram. However, 
Patel et aP find it to be in accord with this prediction, 
after accounting for "experimental noise" which is de- 
termined by measuring the magnetic field asymmetry of 
the CB spacings^ (which indicates motion of the dopant 
atoms, or some other scrambling of the single-particle 
potential). Thus, the Patel et al d&idP seem to lie in the 
weak-coupling region. Accounting for experimental noise 
could also put the Sivan et al data in the weak-coupling 
region, consistent with the equality of parameters in the 
two experiments. 

The experiments by Simmel et al and Abusch-Magder 
et aP are performed on Si quantum dots. Including val- 
ley degeneracy one finds « 2.2, and g « 18. While one 
cannot directly relate rg to Um, one expects that some 
Um might become more negative as increases. (For 
Ts > 2 local charge density correlations develop in the 
dot^**, similar to the classical limit^^. While a Fermi sur- 
face distortion is not a charge density wave, it enhances 
the susceptibility for one, and could thus be a precur- 
sor). Indeed, two signatures of the critical fan are found 
in this experiment'^. The CB peak-spacing distribution 
is found to be 7-8 times wider than expected from the 
Universal Hamiltonian (assuming spin degeneracy). The 
most striking feature of the data is that the width of the 
CB peaks does not vanish'^ as T ^ 0. This is just what 
is expected for a system was located in the critical fan: 
The ground and low-lying excited states are "Fock-space 
delocalized" , and single-particle states are broad even at 
low energy. 

There are many problems one can attack using our ap- 
proach. The d = 2 disordered, interacting bulk system, 
which shows experimental^^ and theoretical^'' indications 
of undergoing a metal-insulator transition, with a Stoner- 
type instability that occurs for arbitrary negative triplet 
Mot^^, can be modeled as a d = 2 lattice of quantum dots. 
In the strong-coupling phase, each dot will have a slow 
collective degree of freedom, which could lead to phase- 
breaking and lack of coherent backscattering of quasipar- 
ticles at very low energies. Persistent current measure- 
ments in systems with a flux, which have resisted the- 
oretical explanation^^, may succumb to our treatment, 
as may disordered gapless superconductors exhibiting a 
novel metal-insulator transition'^^ , since the eigenvalue 
and eigenvector statistics of their grains are expected to 
be governed by one of four newly discovered RMT uni- 
versality classes'^*. Effects of finite T on the CB peak 
spacing'^^ in the strong-coupling phase constitute a nat- 
ural extension of our work. 
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In summary, by identifying a small parameter 1/g, 
we have managed to control the problem of electrons 
subject to interactions and disorder in ballistic quan- 
tum dots. At 1/g = the (exact) saddle-point result 
shows that order parameter cr acquires a nonzero aver- 
age past a critical coupling, and the correlation length 
exponent is i/ = 1, results coinciding with MM's results 
based on the one-loop /3-function^^. The transition is 
the disordered version of the Pomeranchuk transition of 
the Fermi surface^^'^^'^^. For 1/g > Q the critical point 
controls a critical fan in which a fluctuates wildly, lead- 
ing to a ground state not originating from a single Slater 
determinant, and a broad quasiparticle peak even at low 
energy (which has been seen in experiment''). Thus Fock- 
space delocalization is just cr delocalization. The nonzero 
(cr) ~ 0{g) in the strong-coupling phase is seen to ex- 
plain the significant dependence of the CB peak spac- 
ing distribution on the interaction, just as (tr) = in 
the weak-coupling regime makes the distribution insen- 
sitive to interaction. Finally, note that the increase of 
the effective level spacing leads to an enhancement of the 
persistent current^^. 
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